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Abstract. We consider random Fibonacci sequences given by x n +i = 
±/3x n + x n -i. Viswanath (|4|), following Furstenberg (|2|) showed that 
when = 1, linin^oo l^nl 1 /™ = 1.13 . . . , but his proof involves the 
use of floating point computer calculations. We give a completely ele- 
mentary proof that 1.23375 > (E(\x n \)) 1 / n > 1.12095 where E(\x n \) 
is the expected value for the absolute value of the nth term in a random 
Fibonacci sequence. We compute this expected value using recurrence 
relations which bound the sum of all possible nth terms for such se- 
quences. In addition, we give upper ands lower bounds for the second 
moment of the |af n |. Finally, we consider the conjecture of Embree and 
Trefethen (| 1 1), derived using computational calculations, that for values 
of P < 0.702585 such sequences decay. We show that as /3 decreases, 
the critical value where growth can change to decay is in fact -^=. 

Consider the random Fibonacci sequence generated by the recursive rule 
x n+ i = ±x n + x n _i. Clearly, if one considers the expected value E(x n ), 
this will be equal to x , but if one considers instead E(\x n \) this is not the 
case. Here, E(\x n \) is the expected value for the absolute value of the nth 
term in a random Fibonacci sequence. While it might seem possible that 
somehow plusses and minuses cancel, and the average random Fibonacci 
sequence remains bounded, Viswanath (H) showed that in fact almost ev- 
ery sequence grows exponentially with mean value w (1.13 . . .)". For more 
general random sequences x n+ i = ±fix n + Embree and Trefethen 

gave numerical evidence for the dividing line between sequences which 
grow and those which decay (HI) - this case has also been investigated by 
Sire and Krapivsky (Q). 

Unfortunately, Viswanath's proof is quite complex, and involves the use 
of floating point computer calculations. In this note, we give a proof that is 
completely elementary, and while we don't obtain as precise a value for the 
growth rate, we do obtain bounds on the variance for well. 

We consider the set of possible random Fibonacci sequences as a tree, 
with the first leaf being x 2 (which we label n = 0), and each leaf Xi having 
two branches leading to the child leaves x^\ + and x^\ — X{. Such a 
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tree lists the 2 n such sequences of length n. At any given level in the tree, 
we can replace negative values with their absolute values which amounts 
to nothing more than a reversal of the orientation of the leaves (which is 
addition and which subtraction). We will develop recursive formulas for 
the sum of the absolute values of the entries in a row in terms of the sums of 
the previous three rows - the reorientation will allow us to assume that the 
values in our initial two rows are positive. The expected value ^(jx™!) can 
then be bounded since it is simply the sum of the nth row divided by 2 n . 

We begin by considering the leaves below a given entry a in the tree (see 
Figured]). a has two children and there are two possibilities for each, either 
b > a or b < a . It is clearly impossible that both children can be less than 
the parent, and we assume that b\ > a, but we make no assumptions yet 
about the value of b 2 • 

a 

bl b2 

/\ /\ 

ci = la + bil di= |a — bil C2 = la + b2l d2 = la-b2l 

= a + bi =bi-a =a + b2 



bl +ci 



ci-bi 



bl + dl bl -dl 



b2 + c2 c2-b2 



b2 + d2 



Ib2 - d2l 
= lb2-la-b2ll 



Figure 1 . The tree below a. We assume b\ > a. 



We consider the sum of the last row, and bound this above and below by 
multiples of the sums of the previous three rows. We note that clearly it 
would be possible to eliminate any c's and d's from Figure [H but bounds 
obtained in this way are not as good as the bounds we give below. The 
only term in the sum of the bottom row which depends on the value of b 2 is 
\b 2 — \a — b 2 \\. 

The sum of the last row is given by 

a = a + bi + b 2 + +2c : + 2c 2 + di + d 2 + \b 2 — \a — b 2 \\. 
We have the following upper and lower bounds for this sum, 
Lemma 0.1. 

4a + 61 + b 2 + ci + c 2 + d x + d 2 < a < 4a + 2b x + 2b 2 + c 1 + c 2 + d 1 + d 2 . 



AN ELEMENTARY PROOF THAT RANDOM FIBONACCI SEQUENCES GROW EXPONENTIALLY 

Proof. The left hand inequality is immediate, since c\ + c 2 = 2a + b\ + b 2 
and we assume that 61 > a. To show the right hand inequality, we need 
only show that \b 2 — \a — b 2 \ | < a. We need to consider three cases. 

(1) b 2 > a. Then \b 2 — \a — b 2 \ \ = \b 2 — (b 2 — a)\ = a. 

(2) a > b 2 > a/2. Then \b 2 — \a — b 2 \ \ = \b 2 — (a — b 2 )\ = \2b 2 — a\ = 
2b 2 — a < a. 

(3) b 2 < a/2. Then \b 2 — \a — b 2 \\ = \b 2 — (a — b 2 )\ = \2b 2 — a\ = 
a — 2b 2 < a. 

□ 

If we denote the sum of the nth row by S[n), then the lower bound in 
Lemma lOTl can be written as 

AS[n - 3] + S[n - 2] + S[n - 1] 

while the upper bound is 

AS[n - 3] + 2S[n - 2] + 3[n - 1]. 

These give recurrence relations for the sums of the rows which will give up- 
per and lower bounds for the actual sums. If we assume that our sequence of 
sums consists of rational numbers, solutions to the relations can be obtained 
which depend on the real roots of the irreducible cubics x 3 — x 2 — x — 4 and 
x 3 — x 2 — 2x — 4. The growth rates for S[n] will be given by these roots 
divided by 2 (to account for the doubling of the number of entries in each 
row. 
We get 

Theorem 0.1. 

1.23375 . . . > (E(\x n \))n > 1.12095 . . . 

To determine bounds for the variance, we need to consider the sum of the 
squares of a given row (the second raw moment fi 2 ), which we label SS[n). 
This is considerably easier. If we consider starting again with some entry 
a, with children bi and b 2 , then the grandchildren of a will have absolute 
values I a — &i|,|a + &i|,|a — b 2 \,\a + b 2 \- Squaring eliminates all the cross 
terms, and we just get a sum of squares of 4a 2 + 26 2 + 262, which is four 
times the sum of the squares of the first row plus twice the sum of squares 
of the second. The recurrence relation 

SS[n] = 2SS[n - 1] + ASS[n - 2] 

gives a growth factor for u 2 of 1 + yh. Now, the variance is given by 
« 2 = —Hi + /4 an d we have upper and lower bounds on /i[ using the 
growth factors in Table 1. In either case, as n — > 00, the contribution from 
the fif disappears. This gives 
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Theorem 0.2. 

lim (/i 2 |a; n |)™ = 1 + V5 

n—*oo 

Finally, if we generalize the allowable sequences by allowing multiplica- 
tion by a factor (3, we get a modified version of Figured shown in Figure 
El Here we consider only half the tree below a, and account for possible 
combinations of these halves below. There are two possible sums for the 



a 

i 
b 



\a-(3b\ 



d = \a + (3b\ 



\b-/3c\ 



\b + (3c\ \b-pd\ 



Figure 2. Half the tree below a when a multiplicative con- 
stant is included. 

absolute values of the bottom row, depending on the sign in the absolute 
value on the left in the previous row, 

(1) j3b > a /3(c + d) + 2b + \b + (3a - (3 2 b\ + \b-(3a- (3 2 b\ 

(2) f3b < a (3(c + d) + 26+ \b + (3 2 b-(3a\ + \b - (3 2 b - f3a\ 

There are now various possibilities for these sums, depending on whether 
the number in either or both absolute values is negative. This gives six 
possible sums, listed with the conditions which generate them and the re- 
strictions they impose in Table UJ 

Conditions Restrictions Sum 



1. (3b>a,f3a + (3 2 b > b (3 2 > \ f3(c + d) + 2b + 2(3a 

2. f3b > a, [3 2 b > b + (3a (3 > 1 (3(c + d) + (2 + 2(3 2 )b 

3. (3b > a, (3a + (3 2 b < b (3<l (3{c + d) + (4 - 2(3 2 b 

4. (3b <a,(3a + (3 2 b >b (3{c + d) + {2 + 2(3 2 )b 

5. (3b < a, (3a > (3 2 b + b b = (3{c + d) + 2(3a 

6. (3b <a,(3a + (3 2 b<b (3 2 < \ (3{c + d) + 2b- 2(3a 

Table 1. 
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Note that if one combines these half sums, a number of combinations are 
impossible because of the conditions for such a sum to occur. It is possible 
that both children of a will satisfy one set of conditions, thus combinations 
do include some use of the same row twice. Because of the subtraction, the 
sixth row must be involved if the expected value is to decay exponentially, 
and this means that exponential decay can only occur if the condition for 
this is satisfied, namely if 

B 2 < - 3 < ^ ~ 0.707107. 

One would not expect decay to begin immediately as 3 crosses this value, 
since the growth or decay rate will be a combination of various possible 
relations. x 

We show graphs for some other possible growth factors for values of 3 
between and 1 in Figure |3] If, instead, one fixes a row index n, it is 
clear that the expected value of the sum of the nth row does not depend 
smoothly on 3, although one needs to "zoom in" quite far to see this for 
any reasonable large values of n (Figure |4|). It seems possible that further 
investigation of the probability for each row to occur might allow for a 
deeper understanding of such scaled random Fibonacci sequences. 

Growth Factor 



Growth Factor 




Rows 3 and 6 form Table 1 . Rows 4 and 6 from Table 1 . 



Figure 3. Possible growth and decay as a function of 3. 
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Figure 4. Mean value of the row sum for a fixed row index n. 
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